A common solution to the cosmic ray anisotropy and gradient problems 
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Multichannel Cosmic Ray (CR) spectra and the large scale CR anisotropy can hardly be made 
compatible in the framework of conventional isotropic and homogeneous propagation models. These 
models also have problems explaining the longitude distribution and the radial emissivity gradient of 
the 7-ray galactic interstellar emission. We argue here that accounting for a well physically motivated 
correlation between the CR escape time and the spatially dependent magnetic turbulence power can 
naturally solve both problems. Indeed, by exploiting this correlation we find propagation models 
that fit a wide set of CR spectra, and consistently reproduce the CR anisotropy in the energy range 
lO'^ — 10* GeV and the 7-ray longitude distribution recently measured by Fermi-LAT. 
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Introduction : The propagation of Cosmic Rays (CR) 
throughout the Galaxy is far from being fully understood. 
One of the most important, and still uncertain, quantities 
which determine the CR spectra and spatial distribution 
in the Galaxy is the power-law spectral index [6) of their 
rate of escape as a function of energy. This quantity is 
typically probed by secondary /primary CR ratios, the 
Boron to Carbon ratio (B/C) most commonly, which is 
measured with high accuracy up ~ 1 TeV. Several de- 
generacies with other physically relevant processes, how- 
ever, prevent from casting stringent constraints on 5. In- 
deed, depending on the adopted values of other (poorly 
known) parameters as the Alfven and convective veloci- 
ties, we can have 5 ^ 0.5—0.6 for plain diffusion (PD) and 
low reacceleration models (see e.g. [l|-|3[), 5 ~ 0.8 — 0.9 
for convective models [J| or 5 ~ 0.3 for non-convective - 
strong reacceleration models |5|. In principle, measure- 
ments taken in other CR channels could help reducing 
some of those degeneracies. However, multi-channel anal- 
ysis opens new problems. For example, strong reacceler- 
ation - low 5 models are disfavored by antiproton data 
(e.g. [SJ, y and Ref.s therein). The observed low energy 
CR electron and positron spectra [7| and the spectrum of 
the synchrotron emission of the Galaxy [81 are also hardly 
compatible with those models. Furthermore models with 
small 5 are in tension with the modeling of CR accelera- 
tion in sources as they require too steep source spectra in 
order to reproduce the observed propagated CR spectra 
[9|. On the other hand, models with 5 > 0.5, which are 
preferred on the basis of CR spectral data, face major 
problems with the observed large scale CR anisotropy. 

Diffusion of charged CRs in the turbulent Galactic 
Magnetic Fields (GMFs) erases almost completely the 
information about source positions, leaving only a small 
residual anisotropy. Current data (see e.g. [lOj and ref.s 
therein) show a dipole-like large scale anisotropy (LSA) 
of the order of 10"'^ very weakly energy dependent above 



'^ 10 TeV, plus several smaller scale anisotropics in vari- 
ous directions. Below few TeV the observed LSA depends 
more strongly on energy and drops down to ^ 10~^ at 
100 GcV. Apart the Compton-Gctting effect (yielding an 
energy independent anisotropy at the ~ 10^^ level), the 
origin of the observed anisotropy can be twofold: the CR 
drift due to the inhomogeneous source distribution and 
the global leakage from the Galaxy, and the stochastic ef- 
fect of local sources. While source stochasticity can likely 
explain the observed behavior of the anisotropy above 10 
TeV [I3I , it cannot reconcile the anisotropy data with the 
LSA predicted by diffusive models with S > 0.5 [l3|. We 
refer to such discrepancy as the CR anisotropy problem. 

Non-local observables, like the galactic 7-ray interstel- 
lar emission, are a precious probe of the large scale CR 
spatial distribution, hence they may shed light on the ori- 
gin of the anisotropy problem. The interpretation of the 
angular distribution of that emission in terms of standard 
diffusion models also faces a long withstand ing problem. 
It has been known since the EGRET era [IJ] that the 
CR galactocentric radial distribution, derived assuming a 
source density deduced from pulsar or supernova remnant 
(SNR) catalogues, is much steeper than the one inferred 
from the 7-ray interstellar emission along the Galactic 
plane. This is known as the CR gradient problem. 

Thanks to its high angular resolution, Fermi-LAT, 
upon recently confirming the existence of the gradient 
problem [15|, [ig, allowed also to separate the emission 
coming from CR interactions with the molecular gas 
(whose modeling is strongly affected by the uncertainty 
on the conversion factor between the tracer CO inten- 
sity and the H2 column density, Xco) from the emis- 
sion due to CR interactions with the atomic gas (whose 
density is better known from its 21 cm radio emission). 
An analysis based on 7 -ray maps of the third Galactic 
quadrant [l6| (see also [17|) showed that the emissivity 
from neutral gas shows a radial decrease smaller than the 



one predicted by conventional models. This confirms the 
gradient problem independently of the uncertainties on 
the Xco parameter, thereby strongly disfavoring expla- 
nations invoking a sharp rise of the Xco in the outer 
Galaxy [l8|. Some alternative explanations explored in 
|16| do not appear satisfactory: a thick halo is disfavored 
both by ^°Be/^Be |6| and radio data (see e.g. [19| and 
ref.s therein); a smooth source distribution is in contrast 
with SNR catalogues. Furthermore, because the LSA is 
proportional to the diffusion coefficient (DC) D (see be- 
low) and D/ H (where H is the diffusive halo half-height) 
is fixed against the B/C ratio, a thick halo would worsen 
the anisotropy problem. 

In this letter we argue that the CR gradient and 
anisotropy problems share the same origin, residing in 
a too simplified treatment of CR diffusion commonly as- 
sumed to be isotropic and spatially homogeneous. We 
drop both these assumptions (supported neither by the- 
oretical arguments nor by numerical simulations) and 
show that, by making CR escape faster in the most ac- 
tive regions, hence the CR radial distribution smoother, 
we solve naturally both the CR gradient and anisotropy 
problems. We exhibit physically motivated propagation 
models that, while fitting local data of CR spectra, yield 
also the radial CR distribution required to reproduce the 
diffuse 7-ray emission measured by Fermi-LAT, similarly 
to what we already found to interpret EGRET data [20[ 
(see also 2l|, [2^ for other possible approaches) . More- 
over we show for the first time that these models pre- 
dict a significantly smaller CR anisotropy than the cor- 
responding models assuming isotropic and uniform diffu- 
sion. As a consequence, the agreement with the measured 
anisotropy is very much improved, especially in the en- 
ergy range 0.1 -^ 10 TeV, hinting at a new self-consistent 
and comprehensive description of present CR data. 
Anisotropic inhomogeneous diffusion: Charged CRs dif- 
fuse in the turbulent component of the GMF. The regular 
GMF (directed almost azimuthally along the spiral arms) 
breaks the isotropy hence the diffusion tensor has to be 
expressed in terms of two coefficients D\\ and D±^ describ- 
ing diffusion in the parallel and perpendicular direction 
to the regular GMF. For weak turbulence rj = 6B/B <C 1 
{SB represents the rms fluctuation of the magnetic field) 
quasi linear theory (QLT) predicts D± « ^\\V^^ hence 
perpendicular diffusion is irrelevant in that case; in the 
opposite case isotropy should be restored. 

In the Galaxy, however, where the regular and the tur- 
bulent components have similar strengths (77 ~ 1) the 
situation is more complex and can be reliably studied 
only by means of numerical simulations. Here we refer 
to the results reported in [23[, where propagation of CR 
protons in simulated realizations of random MFs with 
r] = 0.5-^2 and in the presence of a large scale azimuthal 
field was discussed. For the random component, a Kol- 
mogorov spectrum was assumed in agreement with obser- 
vations |24| . Simulations were performed only for proton 



energies above 10^^ cV (10^^ eV for ry = 1) due to com- 
puter time limitations. Since the simulated values of -Dm 
and D± decrease slowly and steadily with decreasing en- 
ergy, and there are no reasons why such behavior should 
change at lower energies, we assume that their ratio can 
reliably be extrapolated down to the energies considered 
in this work. As we will show, this assumption is indeed 
consistent with several experimental facts. Remarkably, 
only the simulated value of the D±/D^\ ratio is relevant 
here since the absolute values of the DC components will 

be fixed against CR data. 

Similarly to QLT predictions, the authors of [22] found 
Du and D± to behave oppositely with respect to the tur- 
bulent power, the former (latter) decreasing (increasing) 
with rj. This can be understood as GMF lines random 
walk being enhanced with increasing turbulence strength. 
Moreover, they found that while D\\ ex E^^^ (in agree- 
ment with QLT) the slope of D± is significantly steeper: 
D± oc £;0.5-H0.6^ Perpendicular diffusion should thus be 
the dominant CR escape channel in the inner Galaxy and 
in the next outer regions [23|. Indeed, for H < i?Gai— ^o, 
with i?Gai — 20 kpc and Rq ~ 8.5 kpc, the paral- 
lel/perpendicular escape time ratio is 
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where i?arm — 7ri?0 is roughly the length CR have to 
travel to escape the Galaxy along the spiral arms, parallel 
to the regular GMF. For 77 - 1 - 2, Dj_/D\i > 10"^ 
was found at 10^ GeV slowly decreasing with decreasing 
energy [23[. Extrapolating this result down to ~ 10 GeV, 
we still find Tl < Ty unless H is considerably larger than 
the preferred value iJ ~ 4 kpc determined on the basis 
of the observed ^°Bc/^Be ratio and radio data. 

The above considerations have two main consequences: 
1) Under the obscrvationally preferred conditions of Kol- 
mogorov turbulence, the CR escape time is expected to 
depend on energy as E^-^^^-^. Noticeably, this is the 
same dependence which is favored by a combined anal- 
ysis of CR nuclei and antiproton spectra [3[. 2) In the 
inner Galaxy at least, Tcsc ^ T± = H'^/6D± should be 
anti-correlated to the turbulent power, hence to the den- 
sity of CR sources (see also [20|), which we assume inject 
turbulence in the ISM. 

Solution of the CR gradient problem: We study the ef- 
fects of our assumption of a spatial correlation between 
the DC and the source density by solving the diffu- 
sion equation with the DRAGON numerical diffusion code 
[25| . which, differently from other numerical and semi- 
analytical programs, is designed to account for a spa- 
tially dependent DC. The code is 2-dimensional {R, z) 
and assumes a purely azimuthal (no arms) structure of 
the regular GMF. Therefore we can only model perpen- 
dicular diffusion and the DC is treated as a (position 
dependent) scalar. Nevertheless, as only the escape time 



is relevant to determine the CR density, we can account 
for parallel diffusion along the spiral arms by using an ef- 
fective DC: D^s{R) = max [Dj_{R), {H/R^,^f D\\{R)\. 
We assume therefore the phenomenological dependence 
D±{R) oc Q{Ry, where r > is a free parameter to be 
fixed against data (simulations do not allow to determine 
T with sufficient accuracy). According to QLT and nu- 
merical simulations we assume I?|| to have an opposite 
dependence on the turbulence strength, hence D^^{R) oc 
Q{R)~'^ . We remark that parallel diffusion has almost no 
effect on the 7-ray angular distribution and the local CR 
anisotropy, as it becomes relevant only in the most exter- 
nal regions of the Galaxy, where the source density (hence 
turbulence injection) is very small. Its presence, however, 
naturally prevents the escape time from taking unphys- 
ical large values at large R. For the source radial dis- 
tribution we adopt Q{R) ex (i?/i?0)i-9exp(-5(-^^)), 
based on pulsar catalogues [26|. Using other, observa- 
tionally determined, distributions would not change our 
main results. Similarly to [3|, |20[ we assume a vertical 
profile Dcs{R,z) = Des{R)exT>{z/H). We also assume 
D ex (w/c)"""' {v is the particle velocity) to reproduce the 
low-energy B/C data as shown in those papers. This does 
not affect the results discussed here. We fix 7f = 4 kpc 
and for each value of r we set the D normalization to 
match the observed B/C and other light nuclei ratios. 
We fix the D rigidity dependence 6 = 0.6 in the rest of 
our Letter. To better highlight the effects of inhomo- 
geneous diffusion we consider here only PD propagation 
setups. Adding moderate reacceleration and radially uni- 
form convection does not change significantly any of our 
results. 

We find a good fit of the B/C for all values of t G [0,1]. 
The best fit D normalization only mildly depends on 
T. Also the computed antiproton and mid-latitude 7- 
ray spectra match observations within errors. We then 
calculate the 7-ray emissivity from the CR spatial distri- 
butions in our models. As clear from Fig. [1] the model 
T = (uniform diffusion) does not reproduce the ob- 
served emissivity profile. We obtain the simulated 7-ray 
angular distribution by performing a line-of-sight inte- 
gration of the product of the emissivity times the gas 
density. For consistency we use the same gas distribution 
271 1 and the same catalogue sources [28[ adopted by the 
Fermi-LAT collaboration. We show in Fig. [2] the longi- 
tude profiles of Galactic 7-ray emission and the residuals 
of the models against data for r = and r = 0.85. The 
model r = is clearly too steep compared to data: it 
overshoots the data in the Galactic center region while it 
undershoots observations by several a in the anti-center 
region. Increasing r yields a much smoother behavior of 
the emissivity as function of R (see [l6| for the possible 
reasons why the emissivity in the II and III quadrants 
do not agree entirely). A good match of Fermi-LAT data 
is achieved for r ~ [0.7-^0.9], with r = 0.85 providing 
an optimal fit and improving the residual distribution. 




FIG. 1. Integrated 7-ray emissivity (number of photons emit- 
ted per gas atom per unit time) constrained by Fermi-LAT 
(orange region [Ig|, grey region [1^) compared with our pre- 
dictions for r = 0, 0.2, 0.5, 0.7, 0.8, 0.9 (from top to bottom). 
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FIG. 2. Predicted longitudinal profile of the 7-ray diffuse flux 
along the Galactic plane compared to Fermi-LAT data [2a ]. 
and residuals. Data are integrated over the latitude interval 
jfe| < 5° and in energy between 1104 and 1442 MeV. Solid 
(blue) line r — 0.85, dashed (red) line r = 0. 



Effect on the CR anisotropy: The CR LSA component in 
the radial direction is related to the CR gradient by 
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which we use to compute the contribution of CR dif- 
fusion to the LSA starting from the CR distribution 
computed in the same PD models as in the previous 
section. Remarkably, with increasing t, hence with a 
smoother CR distribution, the predicted LSA also de- 
creases. Changing from r = to t = 1 reduces the 
anisotropy by almost a factor of 10. Intriguingly, we can 
reproduce the CR anisotropy data [l0| up to few TeV 
with r = 0.85. The discrepancy between our model and 




lems in terms of a s pati ally varying convectivc velocity 
was proposed in 2l|, \2^ . A possible consistent solution 



FIG. 3. The CR anisotropy measured by several experiments 
is compared with our predictions for r = (dashed red line) 
and T — 0.85 (solid blue line). Triangle/circle data were taken 
from muon/EAS detectors as reported in [ly| and [ill ]. 



the observed anisotropy above that energy is probably 
due to source stochasticity which we did not account for 
in this work. Indeed, while below 10 TeV the observed 
anisotropy phase (see [10[ and ref.s therein) keeps almost 
constant to a value compatible with expectations from 
the global CR leakage, above that energy it significantly 
fluctuates, as expected if the contribution of stochastic 
sources becomes dominant. 

Conclusions : In this letter we presented a consistent so- 
lution to the CR gradient and anisotropy problems. Our 
approach is based on the physically motivated hypothesis 
that the CR diffusion coefficient is spatially correlated to 
the source density: regions in which star, hence SNR, for- 
mation is stronger are expected to show a stronger turbu- 
lence level and therefore a larger value of the perpendicu- 
lar DC (oppositely to what happens for Z?||). The escape 
of CRs from most active regions is therefore faster, hence 
smoothing out their density through the Galaxy. Corre- 
spondingly, the predicted CR gradient and anisotropy 
are reduced. We implemented a phenomenological real- 
ization of this scenario and checked that - while CR data 
are still correctly reproduced - our approach also gives a 
remarkably good description of the spectrum and longi- 
tude distribution of the diffuse 7-ray emission measured 
by the Fermi-LAT collaboration. Our analysis provides 
for the first time a unified propagation model which re- 
produces local nuclear spectra and also explains non-local 
obscrvables, and in particular reconciles the preferred 
low-reacceleration models with (5 ~ 0.5 hinted at by the 
combined spectra of nuclei (B/C), antiprotons, electrons 
and radio data (and phenomenologically preferred by ac- 
celeration theory) with anisotropy and gradient obser- 
vations. We take these results as an encouragement to 
pursue a self-consistent theory/computation of non-linear 
CR - MHD turbulence interaction in the Galaxy. We no- 
tice that an alternative solution of the CR gradient prob- 



of CR isotropy problem also deserves to be investigated. 
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